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Abstract 

Extensive magnetization and magnetostriction measurements were carried out on a single crystal of 
GdNi2B2C along the main tetragonal axes. Within the paramagnetic phase, the magnetic and strain 
susceptibilities revealed a weak anisotropy in the exchange couplings and two-ion tetragonal-preserving 
a-strain modes. Within the ordered phase, magnetization and magnetostriction revealed a relatively strong 
orthorhombic distortion mode and rich field-temperature phase diagrams. For H / /{IQQ) phase diagram, 
three field-induced transformations were observed, namely, at: Hd{T), related to the domain alignment; 
Hair), associated with reorientation of the moment towards the c-axis; and Hs{T), defining the saturation 
process wherein the exchange field is completely counterbalanced. On the other hand. For H/ /(OOl) phase 
diagram, only two field-induced transformations were observed, namely at: Hj^iT) and Hs{T). For both 

phase diagrams, Hs{T) follows the relation Hs I - [T/T^f ' kOe with Hs{T 0)==128.5(5) kOe and 

Tn{H = 0)=19.5 K. In contrast, the thermal evolution of Hb{T) along the c-axis (much simpler than along 

the a-axis) follows the relation Hr [1 - T/Tr]^ kOe where Hr{T 0)=33.5(5) kOe and Tr{H = 0)=13.5 
K. It is emphasized that the magnetoelastic interaction and the anisotropic exchange coupling are important 
perturbations and therefore should be explicitly considered if a complete analysis of the magnetic properties 
of the borocarbides is desired. 

PACS numbers: 75.30.Kz,75.30.Cr,75.80.+q 



I. INTRODUCTION 

The field-temperature {H — T) magnetic phase diagrams of the intermetallic borocarbides 
i?Ni2B2C (i?=magnetic rare earth) manifest wide varieties of zero-field magnetic structures such 
as Neel-type, helimagnetic, amplitude-modulated, and squared-up states.^"* Most of these states, 
in particular the modulated ones, are unstable against field and temperature variation, leading to 
a cascade of transformation a^i?i^ similar to the ones observed in the elemental rare-earths. It is the 
general opinion that the remarkable features observed in these H — T phase diagrams are governed 
by the combined influence of exchange couplings and crystalline electric field (CEF) interactions^: 
such an approach reproduced successfully the gross features of the magnetic H — T phase dia- 
grams of HoNi2B2C. Nevertheless, various experimental observations suggest that two additional 
interactions, namely, anisotropic exchange interaction and magnetoelastic coupling, are necessary 
ingredients for the understanding of the behavior of the magnetic borocarbides. Magnetoelastic 
interactions are spectacularly manifested in the onset of a tetragonal-to-orthorhombic distortion 
at Tat for R= Er, Ho, Dy, Tb compoundsiSjiLiSii^, while the anisotropic exchange interactions 
were reported to be necessary for the description of the low-temperature magnetic properties of 
ErNi2B2C (Refi^^) and GdNi2B2C (Ref.^-). Obviously, the anisotropic exchange interactions and 
(two-ion) magnetoelastic couplings can be conveniently investigated in GdNi2B2C wherein CEF 
anisotropy is negligible and the de Gennes factor is the strongest. 

This work reports a study of the magnetic and magnetoelastic properties of single crystal 
GdNi2B2C. The H — T magnetic phase diagrams for fields along the a- and c-axis were determined. 
Moreover, within the paramagnetic phase, the parastriction and paramagnetic measurements re- 
vealed the presence of weakly anisotropic exchange interaction as well as two-ion magnetoelastic 
couplings while within the ordered state these perturbations were observed to be strong enough to 
induce a noticeable modification in the nuclear and magnetic structures. 

GdNi2B2C crystallizes in a body-centered tetragonal structure (space group I4/mmm, point 
symmetry of the Gd site is 4/mmm).^^ The zero-field magnetic state immediately below Tjy = 19.5 
K is a transverse, incommensurate, sine-modulated structure (moments along b-axis) with a wave 
vector ka that decreases linearly from 0.551a* at T/v to 0.550a* at Tr ~13.5 Y^^A^i^^i^^i^ At Tr, a 
moment reorientation sets-in leading to an additional modulated component transversely polarized 
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along the c-axisi^*i^ On further temperature decrease, ka increases monotonicahy reaching 0.553a* 
at 3.5 K. 



II. THEORETICAL BACKGROUNDS 

The ''/^S'-character of GdNi2B2C leads to a greater simplification of its magnetic Hamiltonian 
which is considered to consist mainly of an exchange, a Zeeman, a two-ion magnetoelastic, and 
an elastic term. Within the mean-field approximation (MFA) the tetragonal-invariant bilinear 
exchange interaction can be written as a sum of two terms>i^ 

nT. = -{g.jf^Bfn'^' <J>J, 

= (2 < J, > J,- < J, > J,- <Jy> Jy) (1) 

TC^^ corresponds to the isotropic bilinear exchange coupling while corresponds to an anisotropic 
bilinear coupling. The values of the isotropic (n'^^) and anisotropic (n"^) effective exchange coef- 
ficients can be determined from the paramagnetic susceptibilities, measured with the field along 
(X||c) and perpendicular (x±c) to the c-axis: 

_ C 

C 

^"'^ ~ T-{n^^ +2n"2)C ^ ^ 

C is the Curie constant. In Gd-based compounds, the two- ion magnetoelasticity is related to 
the modification of the magnetic interactions by the strains. For strains that preserve the initial 
tetragonal symmetry, the two-ion magnetoelastic Hamiltonians can be written as>^ 

K}E = -{D^le'^^+D2le'^-)<J>J, 

'^ME = - [Kl^"' + ^Se"') (2 < J. > Jz- <Jx> Jx- <Jy> Jy) (3) 

where D^j are the two-ion magnetoelastic constants and e"' are the normalized, symmetrized 
strains, e"^ is the volume strain while e"^ is the axial strain acting on the c/a ratio: 

e ^ = —={exx + ^yy + ^zz) 

The minimization of the free energy with respect to each strain yields the equilibrium values 
of, say, the ct-strains: 



where M^j (i,j=l,2) is a combinations of the D'^j and the symmetrized elastic constant 
Cq * (/i;=l,2,12). The magnetoelastic couplings can be experimentally probed by measuring the 
macroscopic length change under an applied magnetic field. In a tetragonal symmetry, the relative 
length change, A = measured along the {Pi, 132, f^s) direction when an external field is along the 
(x,y,z) direction is given by: 

^^^-^■^Xxyz = ^e"^ + ^6°^ (2/3| - - Pi) + -^e-TiPl - Pi) (6) 
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where corresponds to the strain mode that leads to a symmetry lowering from tetragonal to 
orthorhombic. This mode can be determined from the appropriate magnetostriction measure- 
ments, namely: 

^7 _ 1 /lOO X 010 A, \ (rj\ 

e — Aioo - Aiooj \J) 

The thermal and field variation of the measurable magnetostrictions arise from the T- and 
i?-dependence of < J,; >. The calculation of such a dependence is usually carried out either by 
proper diagonalization of the total Hamiltonian (mostly appropriate for the ordered state and 
higher applied fields) or by perturbation methods (appropriate within the paramagnetic phase and 
for small magnetic fields). Here we will be concerned only with the paramagnetic case (for more 
details see Refj^^). Utilizing the susceptibility formalisnJ^, one gets the following useful relations: 

^"-^ = X^oiH' (8) 

where (^=li2) is the strain field susceptibility along the indicated direction. A complete 
description of the parastriction can be obtained by combining these equations and Eqs Altogether 
with the relation: 

Xi 

< Ji >= — H; i = x,y, z; Na is the Avogadro number (9) 

NA9^J'B 

III. EXPERIMENTAL 

Single crystals were grown by floating zone method»22 X-ray and magnetic measurements re- 
vealed a single phase of an excellent quality. The size and orientation of the rectangular spark-cut 
bars were chosen to suit the requirement of the magnetization or the magnetostriction measure- 
ments. The isofield and isothermal magnetization curves were measured by the extraction method 
for temperatures down to 1.5 K and fields up to 160 kOe. The sensitivity of the measurement is 
better than 5xl0~^ emu. 

The magnetostriction measurement were performed using a high-accuracy capacitance dilatome- 
ter which allows for measuring thermal expansion and forced magnetostriction in the temperature 
range of 3 - 250 K and under magnetic field up to 65 kOe. The typical resolution is better than 1 

o 

A. In addition, the set-up allows for the rotation of the capacitance cell around the vertical axis of 
the cryostat while the field is maintained in a horizontal direction. In this way, the angle between 
the sample axis and the field direction can be varied across the range — 180°. The relative length 
variation for temperature scan at fixed fields is defined as \h = \f'{T) — 1{Tq)] /1{Tq) while that 
for field scan at fixed temperatures is defined as \t = [liH) — 1{H = 0)] /1{H = 0). Diamagnetic 
contributions and demagnetizing fields are expected to be small and as such no correction for their 
infiuences on the magnetization or magnetostriction curves were attempted. 

IV. RESULTS AND ANALYSIS 

Below we show representative magnetization and magnetostriction curves taken within the 
paramagnetic and the ordered phase. By applying the field along the (001) or (100) direction, 
it was possible to investigate separately the volume distortion e"^, the axial distortion e"^ or the 
orthorhombic distortion . The paramagnetic and parastriction curves were analyzed according 
to susceptibility formalism given in Sec. II. Based on the overall analysis of the ordered state, H — T 
phase diagrams for H/ / (100) and H/ / (001) are constructed. 
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FIG. 1: Thermal evolution of the paramagnetic susceptibilities (circles) and their reciprocals (squares). The 
open and closed symbols represent, respectively, measurements along the a-axis and the c-axis. Each point 
was obtained from an Arrott plot (M^ versus H/M). The effective moment and Curie- Weiss temperature 
for each orientation are given in the inset table. 



A. Paramagnetic Phase 

1. Magnetization 

The paramagnetic susceptibilities were deduced at each temperature (25 K< T <300 K) from 
the Arrott plot (M^ versus H/M). The thermal evolution of these susceptibilities along each of 
the (100) and (001) direction (shown in Fig.l) reveals a weak anisotropic character, similar to the 
one reported by Canfield et ali^. From the analysis of Fig.l and according to EqO we obtained 
n"^=0.12(ib0.04) mole/emu and n"^=-0.05(ib0.04) mole/emu. The experimental error involved in 
the determination of these and n"^ were dictated by the experimental resolution, nonetheless, 
the determination of their signs and their magnitudes are consistent with the results of Canfield 
et alM: on the one hand, n"^ is positive and weak indicative of a small effective 9. On the other 
hand, the magnitude of n"^ is twice smaller (indicative of a weak paramagnetic anisotropy) and 
its sign is negative (i.e. x\\c < X±c) implying that it is energetically favorable for the spins to lie 
within the basal plane. 

The effective moments deduced from the slope of the inverse susceptibilities were found to be 
^e//[ioo]=8.33(4) and /ie//(ooi)=8.12(4)^B; both are slightly higher than the theoretical value 
(;Ue//=7.94/iB). Such a small discrepancy arises possibly due to matrix contribution: an exchange- 
enhancement of the Ni-sublattice susceptibility and/or a polarization of the conduction band. 
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2. Parastriction 



Representative curves of the parastrictions Aioo, Aioo, and Aioo are shown in Fig. 2. 
Similar curves (not shown) were obtained as well for *Aoio and *Aooi {i = (100), (010), (001)). In all 
cases, A is of the order of few 10~^, explaining the earlier report that magnetoelastic effects were 
apparently absent in the temperature-dependent XRD measurementsi^ 
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FIG. 2: Representative parastriction isotherms with the field along the a-axis. The relative length variation 
A is measured along (a) the a-axis, (b) the b-axis, and (c) the c-axis. 



Fig. 3 shows the symmetrized, normalized strains e"^ and e"^ which were constructed from 
parastrictions curves (as those in Fig.2) according to EqEI Both e"^ and e"^ are nonzero and follow 
faithfully the quadratic field dependence (see Eq|Hl). These results coupled with the tetragonal 
character of GdNi2B2C predict a magnetic contribution to the thermal expansion, in addition to 
the conventional contribution expected in the diamagnetic isomorph i?Ni2B2C (ii=Lu,Y, La). The 
symmetry lowering mode e'^ (not shown) is mostly field-independent ruling out any possibility for 
a field-induced orthorhombic distortion within the paramagnetic phase. 

Fig.4 shows the thermal evolution of Xioo 

(T) and 

Xooi(-^)- These curves were obtained from the 
analysis of Fig.3 according to Eq|HJ Following the same procedure, we obtained the x"oo(-^) 
XoQi(r) curves (shown in Fig. 5). Within this paramagnetic range and according to EqsElwe get:- 

"hZ/ioo ^h//ooi 

then: 

ai ai 

AIOO ~ AOOl 

XiSo = XoSi (11) 
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FIG. 3: The normalized symmetrized strains (see Eq0J) plotted against the square of the field {H//a). (a) 
e"^ vs and (b) e"^ vs H^. The continuous lines are the least-square fit of EqlHlto the experimental points 
(symbols). 



Within the experimental errors, these equalities were experimentally confirmed in Figs. 4-5: the 
analysis of the x°'^ curves gave: 

l/^xSom = 361(30) +41(l)xT 

l/^Xooim = 324(73) + 36(2)xr (12) 

while for the curves gave: 

y^/xi^oiT) = -27(17) + 20(l)xr 

l/V^ooim = -37(26) + 21(l)xr (13) 

Thus, the thermal evolutions of the two ai- curves (as well the a2-curves) are similar: the small 
numerical discrepancy is attributed to artefacts stemming from the very small values of x"' = 
1,2). 

From these expressions and the above-cited equations, M^j (i, j=l,2) can be straightforwardly 
calculated. However since the values of the symmetrized elastic constant Cq ''(^=1,2,12) are un- 
available, these calculation are not sufficient for determining the extremely useful magnetoelastic 
constants -D"'-. 
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FIG. 4: The thermal variation of the strain susceptibilities (a) Xioo(^)i Xooi(^) ^^'^ (b) Wxioo(^) ^^^^ 
•\/l/Xooi (2^) • The continuous lines are least-square fitting of the experimental points (symbols) . 




FIG. 5: The thermal variation of the strain susceptibilities (a) Xioo(^)' XooiC^) (^) Vl/XiooC^) ^^'^ 
•\/l/xoQx (7") • The continuous lines are least-square fitting of the experimental points (symbols) . 
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B. The Ordered Phase 

1. Magnetization 

Fig. 6a shows representative magnetization isotherms along the a-axis, Mioo- Each Miqq{H) 
isotherm increases monotonicahy with the field and, for a certain field range, such an increase 
is linear. Moreover, for low temperatures, Mioo(-f^) curves demonstrate a weak event at H£) 
«12 kOe (defined as the field value at which [dM / dH)rp is a maximum) and, on a further field 
increase, a saturation at a relatively high field, for instance //^(l.S K) = 128 kOe. It was observed 
that on increasing the temperature, the saturation process becomes less pronounced. The thermal 
evolution of each of H£) and Hs were plotted in Fig. 11. The saturated moment {7.2^b at T =1.5K) 
is weakly enhanced in comparison with the theoretically expected value {7^b)- This enhancement 
is attributed to matrix contribution, just as the above-mentioned enhancement of the effective 
moment. 




FIG. 6: (a) Representative low-temperature magnetization isotherms along the a-axis (Mioo)- The inset 
shows an expanded view of the low-field part of the magnetization isotherm, (b) The thermal variation of 
the magnetization for representative magnetic field {H//a). The continuous lines in Fig. 6b represent the 
thermal evolution of Hd, Hji, and Hs as determined from the analysis of Figs.6a,8-10b. 



The thermal variation of the isofield Mioo(T) is shown in Fig. 6b. Characteristic features signal- 
ing magnetic events at Hd and Hs were also observed and their values were found to agree with 
the values obtained from the magnetization isotherms of Fig. 6a. This agreement was convincingly 
demonstrated by plotting in Fig. 6b the thermal evolution of each of H£) and Hg- Moreover, Fig. 6b 
shows an extra curve, -ff/j, representing the spin reorientation process (see below). 

Representative magnetization isotherms along the c-axis (Mqoi) are shown in Fig. 7. These 
isotherms are very similar to the Mioo ones except that (i) there are only two events namely at Hn 
and Hs, (h) the field dependence of the magnetization at these events is more pronounced, and (iii) 
Hr was observed to have a different temperature dependence than the corresponding one along 
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the a-axis. However, for both orientation, the saturation is reached at the same Hs{T) curve. The 
thermal evolution of Hs and Hfi deduced from Fig. 7 were plotted in Fig. 12. It is worth noting that 
within the experimental conditions of the extraction technique, no hysteresis features or relaxation 
effects were observed. 




H (kOe) 



FIG. 7: (a) Representative low-temperature magnetization isotherms along the c-axis (A/qoi)- The inset 
shows an expanded view of the low- field part of selected magnetization isotherms, (b) Field derivatives of 
magnetization isotherms showing the onset of Hfi{T) (the field at which the maximum is attained) and 
Hs(T) (the field of maximum inclination). 



2. Magnetostriction 

Fig. 8a shows the forced magnetostriction curves, ^^^Xioo{H), at various temperatures. For 
lower temperatures, a field increase is accompanied by a steep increase in ^'^''Aiooi then a shallow 
maximum, and afterwards a steady decrease with a weak slope. The initial steep increase is very 
likely to be due to the domain-purifying influence of the applied field: the field derivative of A 
shows a maximum at Hd (compare Fig. 6a with Fig. 8b). Surprisingly the ^'^^Aioo curves present a 
very small hysteresis suggesting that the domain distribution has a weak dependence on the sample 
history. 

The thermal evolution of the isofield Xioo{T) (shown in Fig. 9) manifests magnetic events 
that are very much similar (though more pronounced) to the ones observed in the magnetization 
(Figs. 6b) and forced magnetostriction curves (FigS.a). The corresponding field and temperature 
values at these events were plotted in Fig.ll. 

It is interesting to note that for intermediate fields, e.g. H = 15 kOe in Fig. 9, on decreasing 
the temperature below T/v, Xioo{T) shows at first a slow increase till 17.4 K where a decrease 
commences and continues till Tr. Afterwards Xioo(T) increases fast but later tends to saturate to 
1.3x10"^ at r = K. In contrast, for lower (higher) fields, Xiqo{T) is monotonically decreasing 
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FIG. 8: (a) Representative forced magnetostriction curves at various temperatures. Measurements were 
taken while increasing and decreasing the field, (b) Field derivative of these magnetostriction isotherms. 
Note that in (b) the shoulder at the left hand side of HoiT) disappears when the temperature is increased 
above Tji{H — 0) suggesting a correlation between this shoulder and the reorientation process below Tii{H). 



(increasing) with decreasing temperature, tending to saturate at T=0 K. For all applied fields, the 
signature of the event at Tji{H) is clearly observable but tends to fade away as higher fields are 
approached. It is also interesting to compare the magnitude of this two-ion effect with the one 
observed in typical single-ion compound such as ErNi2B2C (Refj^S): for zero field, the saturation 

value of ^°°Aioo in GdNiaBaC is ^-0.5x10"^ while in ErNisBaC is 3x10"^. 

Fig. 10a shows an impressive demonstration of the magnetic domain effects within the ordered 
state. In these curves, the angular dependence of various ^^^Xe isotherms was monitored by rotating 
the single crystal in a field of, say, 30 kOe. Furthermore, Fig. 10b shows the orthorhombic strain e"' 
for H =7.5, 30 kOe, constructed according to EqI7| Under H=30 kOe and a decreasing temperature 
below T/v, e'^ commences to increase slowly till Tr(=15.5 K) and afterwards with a faster rate till 
reaching saturation of -2.6x10" ^. Across the borocarbides, this two-ion orthorhombic distortion 
should vary roughly with the de Gennes factor and as such for ErNiaBaC this partial contribution 
amounts to ~ — 0.4xl0~^ which is one third of the total « -1.4x10"^ observed for ErNiaBaC 
at 2.5 K (Ref.^^): thus for borocarbides, the two-ion magnetoelastic contribution (in particular e'^) 
is not negligible as compared to that of the single- ion effect. 



3. The H ^ T Phase Diagram 

Based on the analysis of the magnetization and magnetostriction curves, we were able to con- 
struct the H — T phase diagrams with the field either along the a-axis (shown in Fig. 11) or 
along the c-axis (shown in Fig. 12). In both diagrams, the boundary describing the saturation 
process (when the exchange field is completely counterbalanced) is well described by the rela- 
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FIG. 9: Representative thermal evolution of the magnetostriction for various fields. Both the field and the 
length-measurement were oriented along the a-axis. The inset shows, in an expanded scale, the thermal 
evolution of the magnetostriction in the neighborhood of Tii{H) and Tn{H). 




FIG. 10: (a) Angular dependence of ^°°X0 isotherms under if = 30 kOe and for T = 3 K, 10 K (below Tr{H)) 
and 16.5 K (above Tr{H)). (b) the thermal evolution of taken at H =7.5 kOe (below Hr{T < 13.5 K)) 
and at H =30 kOe (above Hr{T < 13.5 K)). 
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tion Hs 1 - {T/T^f ' kOe with Hs{T 0)=128.5(5) kOe and Tn{H 0)=19.5 K. The 
second boundary, denoted by Hji{T), is attributed to the reorientation process: its evolution 
for H//c (much simpler than that for H//a) is given by the relation Hji(l-T/TFt)3 kOe where 
Hji{T — > 0)=33.5(5) kOe and Tji{H 0)=13.5 K. The extrapolated values of the characteristic 
points Hs, T/v, are in excellent agreement with the reported onesi ^i^'^i^^i^7i^? The third boundary 
H£)(T), appearing only in the a-axis phase diagram, (see Fig.ll) is related to the domain effects as 
can be inferred from the angular dependence of Fig. 10. a. This H£){T) has the following features: 
(i) the extrapolated Hd{T 0)=12.0(5) kOe while Td{H ^ 0) = Tn{H 0) =19.5 K, (ii) it 
is almost constant for T < 8 K, increases slowly up to a maximum around 16 K and afterwards 
decreases steadily till vanishes at T^: such an increase of Ho with temperature is not usual for 
domain effects. 




T(K) 

Fig. 11 



FIG. 11: The H — T phase diagram along the a-axis as compiled from the magnetization (-^) and mag- 
netostriction (▼) measurements. The continuous line represents the best fit of the experimental points to 
Hs{l-{TlTNf )^ kOe where Hs and Tat were found, respectively, to be 128.5(5) kOe and 19.5 K. 



V. DISCUSSION 

The modulated state in GdNi2B2C (and borocarbides in general) is usually related to the nesting 
features in their electronics structures. Below, we discuss the stability of this state against field 
and temperature variation. First MFA arguments are applied to account for Hs, H^, and the 
linearity of the magnetization isotherms. Next, anisotropic perturbation are introduced and and 
afterwards the anisotropic exchange interaction and magnetoelastic couplings are considered. 

Among the borocarbides, GdNi2B2C has the strongest exchange interactions and the weakest 
magnetic anisotropy. The strength of the exchange interactions (denoted as Hex, or Jeff) is re- 
flected as a high value of Tj^ and Hs (see Figs. 11-12) while the weakness of the magnetic anisotropy 
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FIG. 12: The H~T phase diagram {H/ / (001)) as compiled from the magnetization (V^) and magnetostriction 

(▼) curves. The continuous hne represent the best fit of the data to Hs 1 — (T/Tn)'^ kOe with Hs{T 

0)=128.5(5) kOe and Tn{H -> 0)=19.5 K. The dashed lines represent the best fit to Hr{1-T/Tr)^ kOe 
where Hr{T 0)=33.5(5) kOe and Tr{H 0) = 13.5 K. 



(denoted as Han) is manifested as low Ho and as a weak difference between the parallel and per- 
pendicular paramagnetic susceptibility. MFA arguments can be used to calculate an effective value 
for each of H^x, Jeff, and Han from the observed values of T^r, Hs, and H^. Let us, at first, 
ignore the anisotropic feature {Han=0). Then, from the relation zJeff = 2S(s+i) ^® ™ effec- 
tive nearest neighbors), zJ^ff was calculated ~1.86 K. Substitution of this value into the relation 
Hs = 2zJeffS{S + l)/gfLBS gave Hs ~125 kOe, a reasonable value for the effective saturation 
field. Next, taking into account the case Han 7^ and using Eqs^^ one is able to calculate an 
effective value for each of H^x and Han and furthermore to account for the observed linearity of 
MiooiHn <H<Hs) 

Hs = {Hex — Han) 
Hd = a/ Han Hex — H^^ 

(Mioo) /Ms = H/Hs, for Hd < H < Hs (14) 

Substituting the extrapolated Hs{T — > K) and Hd{T ^ K) (see Fig. 11 below) one obtains the 
satisfactory value Hex{T ^ K) = 129.6(5) kOe and Han{T K)=l.l(3) kOe: as expected for 
a Gd-based compound, the former is much stronger than the latter. 

In the above treatment, the zero-temperature magnetic structure is implicitly assumed as a 
collinear 3d AF state. Rather, for T]i{H = 0) < T < Tjy, the magnetic structure is a modulated 
state with no c-component and the ordered moments are distributed among the a- or b-polarized 
domains. Then applying a field along the a-axis would force the a-polarized domains to align 
perpendicular to the field at H£,{T > T/j) and later on to achieve saturation at Hs{T > Tr). For 
the same temperature range, a field applied along the c-axis is already perpendicular to the domains 
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and as such only the saturation transition would appear. On the other hand for T < Tfi{H = 0), 
the spontaneous emergence of the c-component at Hr{T) constitutes an extra event for both field 
orientations: thus for H/ /a-axis, there are three transitions (at Hd, Hj^, and Hs) while for the 
H/ /c-axis, there are only two transitions (at Hji and Hs)- 

Within the ordered phase, the manifestation of the magnetoelastic coupling and anisotropic 
exchange interactions is much more pronounced, leading to a substantial perturbation of the nu- 
clear and the magnetic structures. To demonstrate the importance of these perturbations, let us 
divide the ordered phase into two regions depending on the magnitude of e'^ . The first region 
is Tr < T < Tjv wherein the low- field e'^ is relatively weak (see Fig. 10) and as such is hardly 
sufficient to modify (substantially) the nuclear structure or the modulated magnetic state. The 
second region is T < Tr wherein the magnitude of the low- field is very large (see Fig. 10) and 
increases whenever the temperature is decreased or the field is increased. As a consequence, the 
nuclear structure is distorted (f < 1 and the four- fold symmetry of the basal plane is reduced to a 
two- fold symmetry). In addition, there is an accompanying dramatic modification of the magnetic 
structure: as the temperature is decreased, the anisotropic character of the exchange interactions 
(just as the magnitude of e'^) is enhanced and the magnetic structure is more and more perturbed 
towards a squared-up state wherein the magnetic moments approach equal amplitudes and pro- 
gressively bunch towards an orientation within the b-c plane. While the square-up character can be 
understood based on entropy arguments, the bunching of the moments is most probably induced 
by the above mentioned perturbations. 

By demonstrating the presence of magnetoelastic and anisotropic exchange couplings in the 
ordered phase of GdNi2B2C, this works confirms the far-sight analysis of Detlefs et al^ and Tomala 
et ali^. On the one hand, Detlefs et al attributed the strong asymmetry in the line shape of the 
resonant XRES below Tr^H = 0) as being due to a straining. Consequently, their XRES spectra 
were analyzed as being composed of two peaks (thus the asymmetric character) below Tji{H = 0) 
and as a single line above Tji(H = 0). It is assuring to note that the thermal evolution of the 
position and intensity of these XRES peaks correlate very well with e"* (compare Fig. 7 of Refi^ with 
our Fig. 10). Tomala et al, on the other hand, analyzed their ^^^Gd Mossbauer spectra assuming a 
transverse (moments along b-axis) sine-modulated state above Tji{H = 0) and a strongly squared- 
up and bunched state below Tfi{H = 0). 

VI. CONCLUSION 

We presented experimental results on single crystal of GdNi2B2C. The parastriction and param- 
agnetic features, well accounted for within the susceptibility formalism, demonstrated the presence 
of weak anisotropic exchange couplings and two-ion e"^ and e"^ strain modes. The strength of the 
latter coefficients indicate that the strain derivative of the coupling constant at the Brillouin zone 
centre have significant amplitudes. The ordered state, on the other hand, is characterized by very 
rich field and temperature phase diagrams that manifest high values of Tjy and Hs (indicative of 
strong exchange interactions), a spin reorientation process at Hji(T), and a domain-purifying field 
Hjj ■ Anisotropic exchange interactions and relatively strong magnetoelastic couplings (particularly 
the orthorhombic distortion mode) were observed, lending support to the reported modification of 
the magnetic structure below Tr. 

It is worth recalling that the modulated state of GdNi2B2C within the temperature range 
Tr < T < Ti\f is very much similar to that of ErNi2B2C within the range 2 K< T < Tjsf, even 
though the dominant aligning forces operating in each isomorphs are very much different. Since the 
exchange interaction and the magnetoelastic couplings are driven by the indirect RKKY couplings 
and that de Gennes scaling holds well among the heavy borocarbides^ then both the anisotropy 
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exchange interactions and magnetoelastic couplings must be present (though with varying strength 
depending on dc Gennes factor) in other borocarbide magnets. As such, these perturbations must 
be explicitly considered when discussing the magnetic properties (in particular the H — T phase 
diagrams) of these magnets. 
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